The impact of massive neutrinos on nonlinear matter power spectrum 
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We present the first attempt to analytically study the nonlinear matter power spectrum for a mixed 
dark matter (cold dark matter plus neutrinos of total mass ~ O.leV) model based on cosmological 
perturbation theory. The suppression in the power spectrum amplitudes due to massive neutrinos 
is, compared to the linear regime, enhanced in the weakly nonlinear regime where standard linear 
theory ceases to be accurate. We demonstrate that, thanks to this enhanced effect and the gain 
in the range of wavenumbers to which the PT prediction is applicable, the use of such a nonlinear 
model may enable a precision of (T(mi/, tot) ~ 0.07eV in constraining the total neutrino mass for the 
planned galaxy redshift survey, a factor of 2 improvement compared to the linear regime. 



Introduction: The relic neutrinos having finite 
masses cause a characteristic suppression in the growth 
of structure formation on scales below the neutrino free- 
streaming scale Exploring this suppression signa- 
ture from large-scale structure probes can be a vital way 
to constrain the neutrino masses In fact the strin- 
gent constraints, m^^tot ^ 0.2-0.6eV, have been derived 
from the galaxy power spectrum and the Lyman-a forest 
power spectrum 0, 0, 01 ■ Planned galaxy surveys such 
as the Wide-Field Multi-Object Spectrograph (WFMOS) 
survey Q will further allow a high-precision measurement 
of the power spectra and therefore continue to improve 
the sensitivity to neutrino masses (e.g., 0]). 

However, most of the previous work on the subject has 
been based on linear perturbation theory for a mixed 
dark matter (MDM) model (see 0] for a review). Even 
at scales as large as ^ lOOMpc relevant for the neu- 
trino free-streaming scale, recent studies based on numer- 
ical techniques or perturbation theory have shown that 
the impact of nonlinear clustering cannot be ignored for 
high-precision future surveys, while these studies focused 
mainly on the nonlinear effect on the baryon acoustic os- 
cillations (BAOs) in the power spectrum jol. \ld^ . Yet, the 
effects of massive neutrinos are ignored in these studies, 
even though the neutrinos with total mass Si 0.06eV, 
implied from the oscillation experiments, cause a ^ 5% 
suppression in the power spectrum amplitude that sur- 
passes the expected measurement accuracy ('^ 1%) at 
each wavenumber band. Also unclear is how the neutrino 
suppression degrades the ability of BAO experiments for 
constraining the nature of dark energy as the neutrino 
effect appears at very similar scales to BAOs. 

In this Letter, we develop a new approach to analyt- 
ically study the nonlinear power spectrum for a MDM 
model, based on perturbation theory (PT). PT is a nat- 
ural extension of the successful linear theory, and is there- 
fore expected to give fairly accurate model predictions up 



to the weakly nonlinear regime. We will then use the PT 
approach to study the impact of massive neutrinos on the 
nonlinear power spectrum, and discuss how the use of the 
PT prediction may help constrain the neutrino masses for 
WFMOS-like surveys. 

Methodology: To develop a PT approach for a MDM 
model, we have to deal with the perturbations of multi- 
fluid components, cold dark matter (CDM), baryon and 
massive neutrinos, which are coupled to each other via 
gravity at redshifts of interest. Hence the expansion pa- 
rameter of PT is not a single quantity in contrast to the 
case of a CDM dominated model in which the expansion 
parameter is only the amplitude of the CDM perturba- 
tions. The density perturbation field of total matter is 
defined as as dm = [Spc + 6ph + Sp^)/pm = fchSch + US^, 
where the subscripts 'm', 'c', 'b', V and 'cb' stand for to- 
tal matter, CDM, baryon, massive neutrinos, and CDM 
plus baryon, respectively, and the coefficients, /cb and 
/j/, are the fractional contributions to the matter den- 
sity, r^ino: fi' = ili/o/^^mo = TOi/,tot/(94.iri,no/i2 eV) and 
/cb = ^ ^ fiy- The total matter power spectrum, Pm(fc), 
is then given by 

Pmik) = /2b^cb(fc) + 2/cb/.Pci,..(fc) + fuPuik), (1) 

where the power spectra with the superscript 'L' denote 
the linear-order spectra (see below) and ^ is the cross 
spectrum between Scb and S^. 

Mixture of the neutrinos in total matter affects the 
nonlinear power spectrum as follows. The neutrinos 
would tend to remain in the linear regime rather than 
going into the nonlinear stage together with CDM and 
baryon, due to the large free-streaming. In addition, the 
prefactor /^ appearing in Eq. ^ is likely to be small for a 
realistic model (e.g. fi, ^ 0.07 in j^]), allowing the nonlin- 
ear neutrino perturbations to be approximately ignored. 
In the following we will thus include only the linear-order 
neutrino perturbations, which can be accurately com- 
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FIG. 1: Fractional difference between tfie mass power spectra 
at 2 = 3 witfi and without the massive neutrino contribu- 
tions, where the two cases fv = O.Of and 0.02 are considered. 
The solid and dotted curves show the PT and linear theory 
results, respectively. The two vertical lines indicate a maxi- 
mum wavenumber limit fcmax up to which the two models are 
expected to be valid (see text). The shaded boxes show the 
expected I-cr errors on the power spectrum measurement for 
the 2 ~ 3 WFMOS survey and the case of f„ = 0.01. 



puted by solving the linearized Boltzmann equations 11 1 . 
The vahdity of our assumption will be shown in [l3 | ^ . 

Following the standard PT approach the CDM 
plus baryon component can be treated as a pressure-less 
and irrotational fluid for the scales of interest. The fluid 
equations for mass and momentum conservation and the 
Poisson equation fully describe the dynamics of the den- 
sity perturbation field, (5cb, and the velocity divergence 
field, 6'cb = V ■ Vch/{aH). The solutions to this sys- 
tem can be obtained by making a perturbative expansion. 
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^c. = sghsT+6T+. ■ . and 9., = eg^+eg^+eHK- . . , 

where the superscript '(i)' denotes the i-ih order pertur- 
bation. In our setting, the nonlinear correction to the 
total matter power spectrum Pm(^) arises only through 
Pcb(^) in Eq- (HI)- The nonlinear Pcb including the next- 
to-leading order corrections is expressed as 
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where the last two terms describe the nonlinear correc- 
tions, the so-called one-loop corrections, that include con- 
tributions up to the third-order perturbations. 



^ In brief we have approximately estimated the nonlinear neutrino 
perturbations by solving the modified Boltzmann equations into 
which the nonlinear gravitational potential including the contri- 
bution of the nonlinear 5cb given by Eq. ^ is inserted, motivated 
by the fact that the nonlinear gravitational clustering is mainly 
driven by the CDM plus baryon perturbations. As a result, the 
neutrino density perturbation is found to be enhanced only by 
up to ~ 10% for fi, Si, 0.05 at scales of interest, corresponding 
to less than 0.01% error in the nonlinear power spectrum ampli- 
tudes due to the additional small prefactor in Eq. JTJ. 



The neutrinos affect the spectrum Pcb through the 
effect on the linear growth rate IJ]. At wavenumbers 



smaller than the neutrino free-streaming scale, fcfs(z) 
0.023(m^/0.1eV)[2/(l + 2)]i/2(f7mo/0.23)i/2 /iMpc"!, 
the neutrinos can cluster together with CDM and 
baryon. Conversely, at fc > kfs , the growth rate of CDM 
perturbations is suppressed due to the weaker gravita- 
tional force caused by the lack of neutrino perturbations. 
Thus the growth rate, Dch{z,k), has a characteristic 
scale-dependence in a MDM model. This fact causes 
one complication in computing the second- and third- 
order solutions for 6ch and 6ch- The /c-dependence of 
Dch causes mode-couplings between the perturbations 
of different wavenumbers in the nonlinear regime in 
addition to the mode-couplings via the transfer function. 
Interestingly, however, we have found that, using the 
analytic fitting formula for Deb in [l3|, this additional 
mode-coupling can be safely ignored for the expected 
small value of fi, [l3|. As a result, the nonlinear spectra, 
^cb^^ ^cb^' ' written in the form similar to that 
for a CDM model case [11: 
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Note that P^^ and P^^^ are roughly proportional to 
the square of P^,^, which enhances the neutrino effect in 
the nonlinear regime, compared to the linear case, P^.^. 

Results; Eqs. ^ and ([3]) show that the PT prediction 
for Pm(fc) at a given redshift can be computed once the 
linear spectra, P^^, P^i „ and Pj", are specified. We use 



the CAMB code [15| to compute the input linear spectra 
for a given MDM model^. Fig. [1] shows the fractional dif- 
ference between the power spectra at redshift z — i with 
and without massive neutrino contributions, where the 
two cases fi, — 0.01 and 0.02 are considered and other pa- 
rameters are fixed to their fiducial values. This plot man- 
ifests several interesting points. First, the massive neu- 
trinos induce a characteristic fc-dependent suppression in 
the spectrum amplitude. For the case of linear theory, 
the suppression becomes nearly independent of k at very 



Our fiducial cosmological parameters are QmO = 0.27 (assuming 
a flat universe), n^-,oh? = 0.1277, r2bo/i^ = 0.0223, = 1, = 
0, = 2.35 X 10~^, and w = —1, where rig, and are 
the primordial power spectrum parameters (tilt, running and the 
normalization parameter) and w is the dark energy equation of 
state. 
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small scales, ^ /cfs, as roughly given by AP/P ^ ~8f^ 
. In contrast the PT result demonstrates that the neu- 
trino suppression is enhanced in the nonlinear regime, 
yielding a new fc-dependence in the spectrum shape. 

Second, comparing the linear theory and PT results 
explicitly tells us the limitation of the linear theory: the 
linear theory is no longer accurate at k ^ 0.2/iMpc~^. 
More precisely, the linear theory result starts to devi- 
ate from the PT result at fe ^ A:L,max = 0.18/iMpc^^ by 
^ 1% in the amplitude, as denoted by the vertical dot- 
ted line^ . However PT also breaks down at scales greater 
than a certain maximum wavenumber limit, /cNL.max, due 
to a stronger mode-coupling arising from the higher-order 
perturbations ignored here. Using iV-body simulations 
for a CDM model, [§] showed that the one-loop PT well 
matches the simulation results up to fcNL.max given by 
A2(fcNL,max,2) = fc^P^ (fc, z)/27r2 1 ^^^^^ ~ 0.4. The 
vertical dot-dashed line denotes fcNL,max derived simply 
assuming this criterion for a MDM model. Thus, in the 
case of z ~ 3, the PT model may allow a factor 4 gain 
in fcmax! observationally, this is roughly equivalent to a 
factor 64(= 4'^) gain in independent Fourier modes of the 
density perturbations probed for a fixed survey volume, 
which in turn improves the precision of the power spec- 
trum measurement. 

Can a future survey be precise enough to measure 
the neutrino effect? This question is partly answered 
in Fig. [H The light-gray shaded boxes around the solid 
curve show the 1-a measurement errors on P(fc) at each 
k bin, expected for the z ~ 3 WFMOS survey (see be- 
low). The neutrino suppression appears to be greater 
than the errors at fc <L 0.03ft.Mpc~^. Another intriguing 
consequence of the nonlinear clustering is that the am- 
plified power of Pm{k) reduces the relative importance of 
the shot noise contamination to the measurement errors. 
This can be seen by the dark-gray shaded boxes showing 
the 1-a errors for the linear spectrum. 

Finally it would be worth noting that wiggles in the 
curves reflect shifts in the BAO peak locations caused 
by the scale-dependent suppression effect due to neutri- 
nos. The amount of the modulations is smaller than the 
measurement errors. Hence the uncertainty in neutrino 
mass is unlikely to largely degrade the power of BAO 
experiments, at least for an expected small [l^ . 

Parameter forecasts: To realize the genuine power 
of future surveys for constraining the neutrino masses, 
we have to carefully take into account parameter degen- 
eracies Q- Here we estimate accuracies of the neutrino 
mass determination using the Fisher matrix formalism. 



^ In Fig. [T] the deviation of dashed curve (linear) and solid curve 
(PT) around fcL.max looks seemingly small due to the fact that, 
for the PT result, the denominator Pj^—q in Pf^/Pj^—Q is also 
computed from PT. 
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FIG. 2: The marginalized 1-a error on the total neutrino mass 
as a function of the maximum wavenumber fcmax used in each 
redshift slice (see text), for the WFMOS survey combined 
with the minimal CMB constraints. The fiducial value of 
= 0.01 is assumed. The solid and dashed curves show the 
results for the PT and linear theory models, respectively. The 
light and dark shaded regions represent the range of k where 
the linear theory and the one-loop PT likely break down due 
to the stronger nonlinearities. 



The observable we consider is the two-dimensional 
galaxy power spectrum given as a function of fc|| and 
k± , the wavenumbers parallel and perpendicular to the 
line-of-sight direction [l6|: 
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where k — (fc^ + fcl)^/^ and ^ = k\\/k. Here, k^_ = 

[DA{z)iid/DA{z)]kf^di_ and fc|| = [H{z)/ H{z)f,d]kf,d\\i 
where Da{z) and H{z) are the comoving angular diam- 
eter distance and Hubble parameter, respectively. The 
quantities with the subscript 'fid' denote the quantities 
estimated assuming a fiducial cosmological model, which 
generally differs from the underlying true model. Al- 
though the equation above simply assumes the linear 
galaxy bias 6i and the linear redshift distortion (3, we 
will instead treat bi and (3 as free parameters in order 
not to derive too optimistic forecasts. This treatment 
would be adequate for our current purpose, which is to 
estimate how PT allows an improvement in the parame- 
ter constraints mainly due to the gain in /cmax. A more 
careful analysis will be presented in detail in [12j . 

Following the Fisher matrix for the galaxy 

power spectrum measurement is computed as 

where pa represents a set of parameters and Cov~^ is 
the inverse of the covariance matrix that depends on 
the power spectrum itself and on survey parameters, 
the comoving survey volume and the number density 
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of galaxies. To compute Fap we need to specify the 
integration range fcmin and fcmax; we will throughout em- 
ploy k^in = 10~^/iMpc~^ to obtain the fully-convergent 
results, and below discuss for the choice of /cmax- Note 
that, for several redshift slices, we simply add the Fisher 
matrices of each slice to obtain the total Fisher matrix. 
The l-cr error on a certain parameter pa marginalized 
over other parameters is given by a^{pa) = (F""^)^^, 
where F^^ is the inverse of Fisher matrix. We employ 
the WFMOS survey parameters in [3| consisting of two 
types of redshift surveys: the 2; ~ 1 survey covering 
0.5 < z < 1.3 with 2000 deg^ and the z 3 survey 
covering 2.5 < z < 3.5 with 300 deg^. We consider 5 
redshift slices. The choice of free parameters is also 
important for the Fisher matrix formalism: we include a 
fairly broad range of the model parameters given by Pa = 
{fJmO, ^^mo'l^ ^^bo'l^ /i., ns, as, A|, , w, f3{zi), bi{zi)}. We 
assume three neutrino species that are totally mass 
degenerate and adopt fi^ = 0.01 as the fiducial value. 
The fiducial ii{zi) and 61 [zj] for the i-th redshift shce 
are computed following [l6|. In total we include 18 free 
parameters. 

Fig. [2] demonstrates the marginalized l-cr errors on the 
total neutrino mass as a function of fcmax, where the 
galaxy power spectrum over a range of fcmin l£ k < fcmax 
is included. The value of fcmax for each redshift slice 
is specified by inverting A^(fcmax;-Zj) for the value given 
in the horizontal axis. The errors shown here are for 
the WFMOS survey combined with the CMB informa- 
tion on cosmological parameters except for the neutrino 
masses, f^, and the dark energy parameter, w. The solid 
and dashed curves show the results for the PT and linear 
theory, respectively. If the linear theory is employed, a 
reliable accuracy to be obtained is cr{m^^tot) — 0.13cV 
in order not to have a biased constraint due to the in- 
accurate model prediction. On the other hand, if the 
PT prediction is vahd up to A^(fcmax) — 0.4 as discussed 
in Fig. [3 the accuracy of <T{m^^tot) — 0.072eV may be 
attainable, a factor of 2 improvement. 

It should be also noted that a wide redshift coverage 
for the planned WFMOS survey is very efficient to break 
parameter degeneracies, especially between the neutrino 
mass and the dark energy parameters 17, [l^, because 
the dark energy is likely to affect gravitational clustering 
only at low redshifts, z ^ 1. 

Discussion: It is of great importance to carefully 
study nonlinear structure formation for a most realistic 
model, i.e. a MDM model including ~ O.leV neutrinos, 
in preparation for future galaxy surveys. While the PT 
model developed in this Letter gives the first step in this 
direction, another complement to the analytic approach 
is to use a hybrid TV-body simulation consisting of cold 
and hot particles, which seems feasible with current nu- 
merical resources by extending the pioneering work [lij 



for a model of ~ lOeV neutrinos to models of ~ O.leV^. 
PT will also play a useful role in calibrating/checking the 
simulations results. 

We have demonstrated that the use of PT may enable 
an improvement in the neutrino mass constraint by a fac- 
tor 2 compared to the case that linear theory is used, for 
the planned WFMOS survey. However our study involves 
several idealizations: most importantly we assumed the 
linear galaxy bias and the linear redshift distortion. At 
least for the large scales ~ lOOMpc, it seems feasible to 
develop a self-consistent model to describe galaxy clus- 
tering observables including the non-linear effects on the 
galaxy bias and redshift distortions for a MDM model, 
by using the perturbation theory approach [2l| and/or 
the halo model approach and by combining with simula- 
tions. Such a refined model to describe galaxy cluster- 
ing observables in the weakly nonlinear regime would be 
worth exploring in order to exploit the full potential of 
the forthcoming galaxy surveys for constraining or even 
determining the neutrino masses. 
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